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ABSTRACT 

Observations by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) are used to quan- 
titatively study the hard X-ray evolution in 5 large solar flares selected for spectral hardening in the 
course of the event. The X-ray bremsstrahlung emission from non-thermal electrons is characterized 
by two spectroscopically distinct phases: impulsive and gradual. The impulsive phase usually con- 
sists of several emission spikes following a soft-hard-soft spectral pattern, whereas the gradual stage 
manifests itself as spectral hardening while the flux slowly decreases. 

Both the soft-hard-soft (impulsive) phase and the hardening (gradual) phase are well described by 
piecewise linear dependence of the photon spectral index on the logarithm of the hard X-ray flux. 
The different linear parts of this relation correspond to different rise and decay phases of emission 
spikes. The temporal evolution of the spectra is compared with the configuration and motion of the 
hard X-ray sources in RHESSI images. 

These observations reveal that the two stages of electron acceleration causing these two different 
behaviors are closely related in space and time. The transition between the impulsive and gradual 
phase is found to be smooth and progressive rather than abrupt. This suggests that they arise because 
of a slow change in a common accelerator rather than being caused by two independent and distinct 
acceleration processes. We propose that the hardening during the decay phase is caused by continuing 
particle acceleration with longer trapping in the accelerator before escape. 
Subject headings: Sun: flares - Sun: X-rays, gamma rays - Acceleration of particles 



1. INTRODUCTION 

Large solar flares are very bright hard X-ray sources. 
The emission originates from energetic electrons with en- 
ergies mainly in the 10s and 100s of keV, believed to 
be accelerated in the corona. These electrons have a 
short lifetime in regions dense enough to generate sub- 
stantial hard X-ray emission, and therefore react quickly 
to changes in the acceleration, transport and emission 
processes. While it may be hard do disentangle the con- 
tributions of the various effects in different flares, the in- 
vestigation of the temporal evolution of the hard X-ray 
spectra in single flares is a valuable tool to study these 
processes. Flare models and theories should be able to 
account for the behavior of the observed hard X-ray spec- 
tra as they change during an event. 

Observations of the spectral hard X-ray evolution have 
revealed two main trends: a soft-hard-soft (SHS) spectral 
evolution of emission peaks, and a progressive hardening 
during whole events (SHH, soft-hard-harder). 

The SHS behavio r of em ission spikes was discovered by 
iParks fc Winckler I (fl969h . and since then has been re- 
ported by many others. Recently. iGrigis fc Benz I (|2004f ) 
surveyed quantitatively the spectral evolution of emis- 
sion spikes during M class events, finding that nearly all 
rise and decay phases of the peaks show the SHS be- 
havior. The excursions in both photon flux measured 
at a fixed energy and spectral hardness can be very dif- 
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ferent from peak to peak. However, they show consis- 
tently a characteristic property: the spectral power-law 
index is a linear funct ion of the logarithm of the flux 
(jGrigis fc Benz Il2005ah . The SHS pattern has also bee n 
observed in looptop sources (|Battaglia fc Benzl [2006*1 . 
Thus it is likely to be a characteristic signature of ac- 
celeration rather than a propagation process. Detailed 
modeling of transit-time damping acceleration of elec- 
trons can reproduce the SHS behavior if the effects of 
particle trapping in t he acceleration regio n and escape 
are taken in account (jGrigis fc Benz H2006h . 

The SHH b ehavio r was first observed by 
iFrost fc Dennis I (|1971[ ). who noted that the spec- 
tral index in the late phase of a flare stayed constant 
at a harder value than measured during the first 
(impulsive) SHS pe ak. F urther even t s wer e studied by 
ICliver et all (|1986l ) and iKiplinger I (|1995| ) using data 
from the Hard X-Ray Burst Spectrometer (HXRBS) on 
SMM. The distinctive feature of the SHH evolution is 
the absence of softening as the flux decreases. 

Kiplingcr found two different subtypes of behavior: 

• hardening during a particular peak 

• hardening during the decay of the whole event 

In the first subtype, substantial hardening occurs dur- 
ing a short period, but after an emission peak the flux 
may soften again. Events of the second subtype (corre- 
sponding to the classic flare with a gradual phase) typ- 
ically have some SHS peaks at the beginning but pro- 
gressively harden afterwards. Despite the name SHH, a 
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general hardening may start already before the largest 
peak. Thus SHH is not limited to the decay phase. We 
note also that the two classes of spectral evolution are 
not clearly separated: most SHH events show impulsive 
SHS peaks in the beginning. The interest in SHH flares 
rose after Kiplinger's report of a high association rate 
between SHH and the occurrence of in terplanetary en- 
ergeti c proton events. More recently, ISaldanha et al. I 
(2008) confirmed the association between solar energetic 
particles and hardenings during the January 2005 solar 
storm events. The link between the two phenomena, 
however, is not the subject of this work. 

The two different kinds of spectral evolution (SHS and 
SHH) seem to support the view that there are two dif- 
ferent stages in flares: an impulsive phase at the be- 
ginning followed by a gradual component, corresponding 
to different acceleration mechanisms. This scenario was 
first p roposed to explain radio observations (|Wild et al. I 
f963). A first (impulsive) phase was suggested to ac- 
celerate electrons producing gyrosynchroton emission, 
and the second (gradual) phase was linked to travel- 
ing shocks (type II radio bursts) accelerating further the 
electrons and also ions. This ide a was then used to in- 
terpret hard X-ray o bservations ([Frost fc Dennis I If 9711 : 
IBai k, Ramatvl[l979l ). Later, the shocks were associated 
with Coronal Mass Ejections (C MEs). Occulted flare s 
seemed to confirm this scenario ([Hudson et al. I fl982^ . 
However, it cannot explain the position of the dominant 
hard X-ray source seen durin g the gradual phas e: imag- 
ing observations by Hinotori ( Ohki et al. I Il983h showed 
that the hard X-ray emission comes from too low in the 
solar coron a to just i fy the con necti o n wit h type II ra- 
dio bursts. iKahler I (|1984f ) and [Ba] ([19861 ) argued that 
the impulsive phase is followed by two independent ac- 
celeration processes. The first happening in the post- 
flare loop arcade is responsible for the late-phase hard 
X-ray emitting electrons, and the second higher up in the 
corona, shock driven, accelerates interplane tary electrons 
and io ns. This was later corroborated by ICliver et al. I 
(1986) using SMM observations. 

Stochastic acceleration reproduces the observed SHS 
behavior but cannot at the same time describe harden- 
ing when the flux decays. If both SHS and SHH phases 
of electron acceleration happen in the same event, why 
does the spectral behavior reverse? Are different acceler- 
ation mechanisms at work, or is there a further parame- 
ter in the same process that changes in the course of the 
flare? To find an observational answer to this question, 
simultaneous imaging and spectral observations are an- 
alyzed that have become available for the first time by 
the Reuven Ram aty High Energy Solar Spectrometric 
Imager (RHESSI: ILin et al. I [20021 ). Its data character- 
ize the spectral evolution of the non-thermal hard X-ray 
flux in unprecedented detail. We also compare the spec- 
tral evolution with the geometrical flare configuration, as 
well as the motions of the coronal X-ray source and the 
non-thermal footpoint sources. 

2. METHOD 

The goal of this paper is a detailed quantitative study 
of the spectral evolution of solar flares showing a hard- 
ening trend in RHESSI observations. Rather than at- 
tempting a statistical study of a large number of flares, 
the analysis is restricted to a few events studied exhaus- 



TABLE 1 

Chronological list of the selected events. 



Date 


GOES peak 


Panel in 




time 


flux 


Figs. [2]and|3] 


7-NOV-2004 


16:05 


X2.0 


F 


10-NOV-2004 


02:13 


X2.6 


A 


17-JAN-2005 


09:52 


X3.9 


B 


19-JAN-2005 


08:23 


X1.4 


C & D 


20-JAN-2005 


07:01 


X7.1 


E 




30 40 50 60 80 100 
Energy (keV) 

Fig. 1. — Example of count spectrum (thick error bars) ob- 
served by RHESSI, integrated for 4 seconds around 17-JAN-2005 
09:43:36 UT. The spectrum was fitted with a thermal component 
(Maxwellian), a non-thermal component (log-parabolic, see Ap- 
pendix A), and the average spectrum before and after the flare, 
assumed to be the background. For clarity, the total model con- 
sisting of the sum of the three components is not plotted over the 
observed count spectra, but its normalized residuals are shown be- 
low instead. The residuals show that the fitted spectral model 
reproduces the observed counts. The reduced \ 2 f° r this spectrum 
equals 0.94. 



tively. Therefore, we do not estimate the occurrence fre- 
quency of SHH flares or their rate of association with 
solar energetic p article events. This has been done by 
Kipli nger 



(|1995h using SMM/HXRBS data, who reports 
24 occurrences of hardening out of 152 events with peak 
flux count rate larger than 5000 counts s . Most of 
the SHH events reported by Kiplinger are in the upper 
M and X GOES class. The reality of this trend needs 
however to be confirmed. 

We selected flares with a GOES flux above XI dur- 
ing RHESSI observation time windows. 50 events satis- 
fying this condition were found in the period from the 
start of the mission (February 2002) to September 2006. 
We additionally required that the rise, main and decay 
phases were well observed to study the spectral evolution 
in time. This left us with 12 candidates. 

As a first approximation, the presence of hardening 
behavior was tested by studying the count rates in the 
energy range from 30 to 60 keV. We fitted the spectral 
index separately in three energy bands (30-40 keV, 40- 
50 keV, 50-60 keV) and looked for either trends of pro- 
gressive hardening or the presence of a late hard phase. 
The lowest of these bands is sometimes contaminated by 
thermal emission. This can be easily spotted by compar- 
ing the time profile with the other bands. Using count 
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Fig. 2. — Temporal evolutions of the spectral index 7, the flux at 50 keV, JgOi and the spectral curvature r} for the selected events. 
Below, the correlation coefficient between the two curves is shown vs. time. It was determined during the time intervals indicated by the 
horizontal bars. The vertical bars represent the 68% confidence range. The onset times of the flare associated CMEs are marked by arrows. 
Ol and 02 indicate linear and quadratic extrapolations, respectively, to the CME altitude vs. time evolution. 

rates is adequate to identify candidate events for spec- 
tral hardening, but may have missed some events where 
hardening happens in a phase of low flux close to the 
background. After discarding two further events with 
high pileup, we found 5 well-observed events with a clear 
signature of hardening. These are listed in Table Q] 

For each of the selected events, the instrumental re- 
sponse matrices and count-rate spectrograms covering 
the energy range from 3 to 500 keV were generated for 
the front segments with a temporal resolution of one 
RHESSI spin period (approximately 4 s). The spatially 
integrated photon spectra in the range 12 to 500 keV 



were fitted with two components: an isothermal com- 
ponent at low energies (below about 20-40 keV) and a 
non-thermal component at higher energies. The shape 
of the non-thermal component is assumed to be a log- 
parabolic curve characterized by three parameters: its 
normalization I Eo at energy eo, its spectral index 7 < 
and its spectral curvature 77. The functional form of this 
model (Eq. IA1[) . as well as the rationale for choosing it 
are explained in detail in App. [A] Figure Q] shows an 
example of an observed count spectrum and the best-fit 
components. The normalized residuals indicate an excel- 
lent fit. 
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10-Nov-2004 



1 7-Jan-2005 





0.001 0.010 0.100 1.000 10.000 100.000 
Flux at 50 keV (photons s 1 cm" 2 keV" 1 ) 

19-Jan-2005 



0.01 0.10 1.00 10.00 100.00 
Flux at 50 keV (photons s 1 cm " 2 keV" 1 ) 



19-Jan-2005 




(D) 



0.01 0.10 1.00 10.00 100.00 
Flux at 50 keV (photons s" cm 3 keV ') 



20-Jan-2005 



0.01 0.10 1.00 

Flux at 50 keV (photons s 1 cm 2 keV _1 ) 

7-NOV-2004 



(E) 




(F) 




0.1 1.0 10.0 100.0 

Flux at 50 keV (photons s" cm" 2 keV' 1 ) 



0.01 0.10 1.00 10.00 100.00 
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Fig. 3. — Evolution of non-thermal X-ray spectral index 7 vs. photon flux /50 in the course of flares. 



During the fitting process, the model photon spectrum 
is folded with the response matrix, yielding the expected 
count spectrum from the model. The background spec- 
trum is then multiplied with a normalization factor A 
and added to the model counts, where A is an additional 
fit parameter for the model fitting and is constrained 
between 0.5 and 2. This correspond approximately to 
the maximum excursion in RHESSFs background in the 
front segments during an orbit. The parameter r\ is 
constrained to be zero or positive (corresponding to a 
parabola in log-log bending down). This ensures that 
the emission approaches for infinitely high energies. 

The large amount of data (more than 3 thousand spec- 
tra) required an automated fitting routine. For every 
spectrum, 2 preliminary passes were done estimating the 
parameters for the thermal and the non-thermal part 
which were then used as starting parameters for the fi- 
nal fitting. This turned out to deliver good fittings for 
most of the data. A check of the quality of the fits was 
performed by looking at the time evolution of x 2 an d of 
the fitting parameters. Spectra with reduced x 2 worse 
than 2 were manually fitted again, and in most cases it 



was possible to find another set of fit parameters yielding 
reduced \ 2 below 2, with the exception of the event of 
20-JAN-2005. 

This event (the largest flare, GOES class X7) is charac- 
terized by very strong thermal emission. At times when 
the non-thermal emission is weak and/or soft, pileup ef- 
fects are especially large in the 20-50 keV band. There- 
fore, the fittings, which are good above 50 keV, have large 
residuals below that energy. This may be due to the fact 
that the photon spectrum model chosen is not suited to 
describe the observed photon spectra, or that the pileup 
correction is inaccurate. Because it is very hard to cor- 
rectly take into account pileup effects in such a regime, 
it is not clear whether the model failure is real or in- 
strumental. Therefore, we let the spectrum model stand 
as it is, but caution that the parameter values fitted in 
the flare of 20-JAN-2005 may be less accurate, due to 
the unknown systematic effects generated by imperfect 
pileup correction. 

For the other events, the fit parameters are of good 
quality and the corresponding photon models are a high- 
fidelity representation of the incoming photon flux. The 
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Fig. 4. — Distribution of th e slope in 7 vs. log(7so) (Fig. [3] 
parameter a in Eqs. Ill and [Bil l for the rise and decay phases in all 
peaks of the observed events (upper panel) and for rise and decay 
phases separately (middle and bottom panel, respectively). 

final distribution of the reduced \ 2 f° r a H events (except 
20-JAN-2007) is such that 89% of all spectra have x 2 
less than 1.5 and 97% of all spectra have % 2 less than 2. 
Therefore the unusual choice of the logarithmic parabolic 
fit-model, explained in Sect. [2 produces good fittings 
and is therefore justified a posteriori. 

Thre e of the events pr e sented here have also been stud- 
ied by ISaldanha et al. I (|2008f ). where the non-thermal 
component has been fitted by a double power-law. The 
temporal evolution of our values for 7 (at 50 keV) is very 
similar to theirs, although the actual numerical values 
are slightly different due to the different fitting methods 
employed and their choice of energy intervals. The dif- 
ferences are larger for the January 20 event, where they 
used the 100-200 keV energy interval. 

The thermal evolution of the events shows a rapid in- 
crease of the emission measure at the beginning of the 
event, followed by a flat peak and a slow decay. The 
temperature is in the range 20-40 MK, peaks before the 
emission measure and decays faster. 

3. SPECTROSCOPY RESULTS 

Figure [2] shows time profiles of the photon spectral in- 
dex 7 and the flux normalization at 50 keV, /50 (given in 
Eq. IA1[) , for the events of Table [IJ as found by the spec- 
tral fitting procedure explained in Sect. [2] and Appendix 
A. 

The observed spectral variability of flares on time 
scales down to ten seconds or less requires the high- 
est possible temporal resolution for the spectral fitting 
(about 4 seconds in our case). Longer integration times, 
while desirable for better photon statistics, are not suit- 
able because the averaging effect of summing spectra 
with different hardness blurs the spectral evolution. Nev- 
ertheless, during the decay phase the flux is so low that 
full resolution spectra deliver noisy values for the fitting 
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parameters. In this case, longer integration times must 
be used. This is acceptable, since the variations of the 
hard X-ray flux are slower during the decay phase, and 
short-lived spikes are less frequent. Therefore, some light 
curves shown in Fig. [2] and subsequent figures use a lower 
cadence of approximately 32 seconds (that is, 8 RHESSI 
rotations) in the decay phase. 

Short gaps lasting about 1 minute can be seen in the 
light curves (Fig. [2]) . They correspond to periods where 
the thick attenuator was removed from the field of view, 
but the X-ray flux was still so large that the dead time 
in the detector prevents meaningful spectral analysis. 

The selected flares show many distinct SHS peaks. 
They are characterized by a temporal correlation of 7 
and log/50, yielding roughly parallel curves in Fig. O 
On the other hand, in the presence of SHH peaks or pro- 
gressive hardenings, the two time profiles diverge. This 
may be illustrated in the event of 19-JAN-2005 (Fig. 
panel C) , where the two profiles run roughly parallel until 
08:29 and then start to diverge. 

To better distinguish between the SHS and the harden- 
ing trends, Fig. [2] also shows the correlation coefficients 
between the spectral index 7 and the logarithm of the 
flux log/50 vs. time. The vertical bars represent the 68% 
confidence range (corresponding to one standard devia- 
tion). SHS peaks are characterized by a correlation co- 
efficient close to +1. Times, on the other hand, when 
the spectrum hardens while the flux becomes lower have 
a negative value of the correlation coefficient. We note 
that during the SHS times correlation is rather constant 
(near +1), whereas during the hardening phase the corre- 
lation coefficients are more erratic and often not close to 
-1. This indicates that it may not be possible to find a be- 
havior similarly well-defined for the periods of hardening 
as it is found for the SHS peaks. In some cases (19-JAN- 
2005, 08:28 to 08:30), the spectral hardness stays nearly 
constant while the flux decays. This yields a correlation 
coefficient near zero. 

It should be noted here that most periods showing 
hardening have low flux. Therefore the corresponding 
time profiles may be more noisy, weakening the correla- 
tion. Some of this effect was compensated by increasing 
the time interval for correlation. Another effect is loss of 
correlation during a broad peak. Again this can be taken 
into account by increasing the correlation interval. 

During the late phase of the event of 19-JAN-2005, in 
the RHESSI orbit following the one featuring the main 
peak, an uninterrupted phase of hardening is seen. From 
09:35 to 10:00 the flux decays exponentially (see Fig. [2l 
panel D). After that time, the emission reaches a hard- 
ness comparable with the one of the background, and it 
becomes impossible to disentangle the two components 
by purely spectral methods. This event will be investi- 
gated in more detail in Section [5] 

We also compared the start of the hardening with 
the onset time of flare-associat ed CMEs, taken from the 
SOHO LASCO CME catalog (jYashiro et al. Il200l . In 
three cases (panel A, B, and F in Fig. [2]), the onset of 
the CME precedes the start of the hardening by 3 to 5 
minutes, in one case by 15-20 minutes (panel C) and in 
one case by 50 minutes (panel E). In all five events an 
associated CME was present, but the hardening phase 
never starts before the CME onset. 

Figure [3] shows the relationship between 7 and log/50. 
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Since the flux increases from left to right and the 
hardness increases from top to bottom, SHS peaks 
show as piecewise linear tre nds with a negative slope 
(jGrigis fc Benz 112004 l2005al ). On the other hand, pro- 
gressive hardening during flux decay is visible as a trend 
with positive slope. Such relations can be written as 

7 = -alog(J eo )+6, (1) 

where a < for the SHS peaks and a > during hard- 
ening phases. The minus sign in front of a takes care of 
the fact that in Fig. [3] the vertical axis for 7 is reversed. 

Careful examination of the plots in Fig. [2] and [3] reveals 
that 

• Most of the emission spikes are well represented by 
straight lines in 7- log /50 during the rise and decay 
phase. The decay phases are sometimes flatter in 
7-log/5o than the corresponding rise phases (e.g. 
panel C), but the opposite is also observed (panel 
A). The event shown in panel E shows some signif- 
icant deviations from the piecewise straight trend. 

• Spectral variability is stronger at the beginning of 
the event (panels A, B, C, F). 

• In the late phase of the events a slower varying 
component is seen, piecewise straight in 7-log/so, 
mostly nearly flat, slowly hardening (panels B, D), 
slowly softening (panels A, F), staying at an ap- 
proximate constant hardness (panel C), or a mix- 
ture of the above (panel E) . 

• During the rise phase up to the strongest peak, 
the hardness tends to increase from peak to peak 
(panels A, C). Events are softer in the beginning 
(all panels). 

Figure [4] shows the distrib ution of the 7-log/5o slopes 
(parameter a in Eqs. [T] and IBljl in the rise and decay 
phases of all emission spikes in the 5 events. If the flux 
increases, a negative slope represents spectral hardening 
and a positive slope a softening. The opposite happens 
during a decay phase. SHS peaks have a < both in the 
rise and decay phase. 

The separate histograms for rise and decay phase dif- 
fer marginally. Spectral hardenings during rise have a 
slightly steeper 7-log /50 slope than softenings during de- 
cays (SHS peaks). Quantitatively, the average value of a 
when restricted to negative values is of —0.70 ±0.06 (the 
uncertainty is the standard error of the mean) during the 
rise phases and —0.51 ±0.08 during the decay phases (af- 
ter removing an outlier with slope -3.2). The average of 
a for the combined set of rise and decays is —0.60 ±0.05, 
corresponding to an average pivot point e nergy (see App. 
B) of e* = 9.4 ± 1.3 in agreement with iGrigis fc Benz I 
l(200l . 

Examination of spectral hardenings during decay and 
softenings during rise phases (trends opposite to SHS), 
restricted to the range 0-3, yield an average a value of 
0.66±0.11 and 0.88±0.18, respectively. These values are 
not significantly different from each other or from the 
corresponding absolute value of the a < averages. 

4. IMAGING RESULTS 

Is there a connection between spectral hardening and 
X-ray source geometry indicating a different acceleration 



site? The position of the hard X-ray footpoint sources 
were investigated in CLEAN images, using detectors 3 to 
8 with a cadence of 60 seconds. In particular, we looked 
for differences in source positions and velocities at the 
onset and during the period of general hardening. 

Figure [5] shows an overview of the position of the ther- 
mal and non-thermal sources during the events. The 
event of 10-NOV-2004 is not shown because it lacks a 
coordinated evolution of source positions. 

The event of 19-JAN-2005 is particularly interesting, 
as it was well observed and the footpoints (FP) clearly 
move alon g the ribbons noticeab le in a TRACE image 
at 1600 A (jSaldanha et al. II2008I ). The motion is fast at 
the beginning and later slows down. The thermal source 
has the form of a loop, rising throughout the event. 

Figure [6] shows the displacement of the northern FP 
source in the directions parallel and perpendicular to the 
northern ribbon 3 . The FP source starts with a velocity 
of 50 km s" 1 along the ribbon and slows down continu- 
ously during the SHS phase until the onset of the hard- 
ening phase, when it becomes nearly stationary. There 
is no evidence of an abrupt transition between the two 
regimes. Furthermore, 19-JAN-2005 is the only event 
clearly showing footpoints drifting apart and coming to 
a stop at the onset of hardening. 

The other events show other positional changes. In 
the following the geometrical behavior of the sources and 
their evolution are described shortly for all the events. 
We distinguish between footpoint sources and coronal 
sources. Because of projection effects, we cannot reli- 
ably reconstruct the three dimensional structure of the 
sources. However we know from limb event observations 
(jBattaglia fc Benzl |2006| ) that FP sources are mainly 
non-thermal and well-observed above 20-30 keV, whereas 
coronal sources are mainly thermal and well-observed be- 
low 20-30 keV. Therefore, we assume in the following 
that the thermal source is coronal and that non-thermal 
sources are footpoints. 

07-NOV-2004: This events features two footpoints 
and a coronal source. The eastern FP moves from W 
to E from 16:20 to 16:24 UT, jumps back to near the 
starting position and moves again from W to E from 
16:25 to 16:30. The western FP moves slightly from SW 
to NE from 16:21 to 16:24, then changes direction with a 
slight jump to W, and slowly moves to NW from 16:25 to 
16:30. The western FP is brighter than the eastern FP 
before 16:24 and dimmer after 16:25. The thermal source 
is located farther E than the FPs and moves slightly to 
N from 16:19 to 16:25, and is not clearly seen in the 
images afterwards, due to the insertion of RHESSI's thick 
attenuators. 

The jump in position around 16:25 roughly coincides 
with the time at which the hardening starts, possibly 
indicating that another loop is actively accelerating 
electrons, but the expansion of the loop, as suggested by 
the FP motions, continues during the hardening phase, 
contrarily to what has been observed in the event of 
19-JAN-2005. 



3 The parallel and perpendicular directions are defined relative to 
the direction of the regression line obtained by least-square fitting 
the positions of the hard X-ray footpoints independently in each 
ribbon. 
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Event Onset of 

hardening 
07-NOV-2004 16:26 

10-NOV-2004 02:12 
17-JAN-2005 09:49 

19- JAN-2005 08:27 

20- .TAN-2005 06:49 

10-NOV-2004: This event has a very complicated 
FP morphology, with sources and source-pairs appearing 
in many different places. It is not possible to find a 
well-defined source motion like in the simpler cases 
with only two footpoints. Here the sources seem to 
jump around as new footpoints in a different position 
become brighter and outshine the old ones. Important 
shifts in position occur at 02:08 and 02:12. The latter 
shift happens at the same time as the onset of hardening. 

17-JAN-2005: Three pairs of FPs are seen. The 
northern pair is stationary and seen from 09:43 to 09:45. 
The southern pair consists of an eastern FP moving to 
SE from 09:45 to 10:05, and a western FP moving to N 
from 09:43 to 09:47, then shifting to W (09:55) and mov- 
ing very slowly to N until 10:05. The last pair of FP is 
to the east and stationary from 10:11 to 10:29. Coronal 
sources are seen in two locations: one to the N of the 
southern FP pair, moving to N from 09:46 to 09:57, and 
the second to the NW of the easternmost FP pair, nearly 
stationary from 10:16 to 10:30. 

There is no clear signature of a discontinuity or a 
change of behavior happening around 09:50, when the 
hardening starts. 

19- JAN-2005: Two footpoints are seen with a 
loop-shaped coronal source between them. The northern 
FP moves to NE from 08:12 to 08:30 (covering nearly 60 
arcseconds) while the southern FP moves, slower, to SE. 
In the meantime, the loop-shaped coronal source moves 
to NW (indicating that it is rising). After 08:30 (only 
3 minutes after the onset of hardening) the northern 
FP is much slower. The coronal source keeps moving 
to NNW. In the next RHESSI orbit (after 09:30), the 
FP sources can still be seen near the old positions at 
08:30. The northern FP is nearly stationary from 09:33 
to 09:59, while the southern FP very slightly moves to 
W, and the coronal source slightly moves to N. 

20- JAN-2005: This near-limb event features two 
FPs and a loop-like coronal source. The southern FP 
moves to W while the northern FP moves to E. The 
eastward motion of the northern FP is not continuous: 
at 06:49 it reverses and recedes until 06:55, when another 
sources appears 20" to E. The double structure lasts 
until 07:01, when the easternmost source fades away. 
The coronal source moves to NW and rises throughout 
the event, slowing down towards the end. The reversal 
coincides with the start of the hardening phase. 

The observed behavior at the onset of hardening for 
all events is reported in Table [2] It summarizes the anal- 
ysis of the source motions in an interval of time span- 



Coronal source 
motion 

uncertain (bad images at 
low energies after onset) 
stationary 
continuous motion 
continuous motion upward 

continuous motion upward 

ning 4 minutes, centered on the onset of hardening. We 
conclude from the imaging observations that there is no 
universal trend holding for all events. Sometimes, there 
seems to be a switch to a different loop system near the 
beginning of the hardening phase. On the other hand, 
such jumps can also be seen during the SHS phase of 
the events, so they need not be significant. There is 
also some indication that the FP motion is slower dur- 
ing the hardening phase, but again this does not hold 
for all events. In the event of 19-JAN-2005, with a sim- 
ple geometry and well observed, the change in spectral 
behavior leading to the hardening phase does not have 
an impact on the morphology of the hard X-ray sources 
seen by RHESSI. 

5. MODELING THE SHH PHASE 

iGrigis &: Benz I (|2006f ) showed that the soft-hard-soft 
trend is expected from a transit-time damping stochastic 
acceleration model that includes escape of particles from 
the accelerator. The hardness is controlled by how fast 
the particle gain energy and how long they are trapped 
in the accelerator. Harder spectra result from longer 
dwelling times of the electrons in the accelerator and 
higher acceleration efficiency. These conditions also al- 
low a larger population of high-energy electrons to build 
up, leading to increased hard X-ray emission from the ac- 
celerator, identified as a part of the looptop source. This 
basic model predicts that harder spectra also have larger 
hard X-ray flux, but cannot explain the soft-hard-harder 
trend seen as the flux decays, because these observations 
associate harder spectra with smaller flux. 

To fit the observed SHS behavior, IGrigis fc Benz I 
(2006]) had to assume that electrons are trapped below a 
certain threshold energy Et and cannot leave the acceler- 
ator. The escaping electron population has a low-energy 
cutoff at Et . Then the photon spectra of the footpoints, 
which dominate the non-thermal emission, harden below 
E^. 

In the following, a simple extension of the basic 
stochastic acceleration model is presented which could 
lead to the observed spectral hardening. We introduce 
the additional assumption that Et increases with time 
in the SHH phase. Therefore, the photon spectrum be- 
low Et hardens while at the same time the flux arriving 
at the footpoints decreases. The important point here is 
that this also happens if the electron spectral index 5 and 
the electron flux normalization in the acceleration region 
are constant. Thus the new variable does not contradict 
the basic properties of the stochastic acceleration model. 
This is shown in Fig. [3 where footpoint photon spectra 
have been computed from a given electron distribution 
with various £/t- As expected, the photon flux decreases 
with increasing Et- Note that the photon spectra have 



TABLE 2 

Source motion near the onset of hardening. 



Footpoints 
motion 

Jump in the position of the eastern FP 

Change in the direction of motion of the western FP 
Jump in the position of the two brighter FPs 
Nearly continuous motion of both FPs 

Continuous motion of the northern FP, slowing after the onset 
Nearly stationary position of the southern FP 

Reversal in the direction of motion of both FPs, slowing down afterwards 
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Fig. 5. — RHESSI hard X-ray positions of the non-thermal footpoint sources (stars) and the thermal coronal source (diamond) for 4 
events. The energy range for the footpoint position is 50-100 keV band (or 25-50 keV at times where the 50-100 keV emission is too faint 
for imaging). The energy range for the coronal source position is 12-25 keV. Systematic motions are indicated by arrows. 



a nearly constant but slightly increasing hardness while 
the flux decreases. 

This extension is compared numerically to the late 
phase of the 19-JAN-2005 event. For simplicity, an elec- 
tron distribution with constant spectral index and flux 
normalization is assumed to escape from the acceleration 
region. It has a low-energy cutoff Et due to trapping in 
the accelerator. Et increases with time. Obviously, we 
do not expect such a simple scenario to reproduce all the 
details of the spectral evolution. The question is, how 
much of the observed features can be explained with this 
simplest extension of the existing acceleration model. 

A non-isotropic energy distribution of fast electrons 
with a positive slope in energy is known to be unstable 
towards growing plasma waves. Therefore our scenario 



also includes an alternative electron distribution featur- 
ing a turnover at Et, that is, a flat distribution below 
Et, instead of a cutoff. The electron distributions with 
cutoff and turnover are assumed as, 



UTOFF 



(E) 



= 1 (£j ifA A| (2) 



if E < E T 



E 



^TURNOVER = ^ / E \ 5 ^ 

F E „ f-^J XE< E T 
The free parameters are the electron spectral index 



Spectral Hardening in Large Solar Flares 



9 




08:20 08:30 08:40 08:50 09:00 

19-Jan-2005 



Fig. 6. — Movement of the northern X-ray footpoint source for 
the event of 19-JAN-2006. The upper curves show the motion in 
the two components parallel (top curve) and perpendicular (middle 
curve, multiplied by a factor of 3) to the northern TRACE 1600 A 
ribbon. The bottom curve displays the cross correlation coefficient 
between hard X-ray spectral index and flux, showing the onset of 
spectral hardening around 08:27. 

r5, the electron flux normalization Fe (electrons s" 1 
keV -1 ), the cutoff or turnover energy Et- The reference 
energy E is fixed at 50 keV. 

An exhaustive description of the method used for the 
comparison of this simple model with the observation 
and finding the best fit is given in App. [Cj 

The late phase of the event of 19-JAN-2005, from 09:32 
to 10:02, comprises 30 minutes of continuous hardening, 
and therefore is well suited for the comparison with the 
variable cutoff model. Figure [5] shows the comparison. 
The observed values of 7 and 77 are plotted as a function 
of /50 together with the best-fit model curve for both 
the cutoff and turnover electron spectra. For the cutoff 
model, the best fit values of the parameters are: spectral 
index S = 6.52, flux normalization at 50 keV F50 = 4.50 • 
10 33 electrons s -1 keV -1 . In the course of the decay, 
Et increases from 98 keV to 159 keV (thus yielding a 
photon spectral index 7 in the range between -2.5 and 
-1.9). The corresponding values for the turnover model 
are: 6 = 6.11, F 50 = 3.67 • 10 33 electrons s" 1 keV" 1 , 
while Et increases from 127 to 224 keV. 

The photon spectra from both the cutoff and turnover 
model are curved downward in the fitted energy range. 
Thus the spectral curvature, 77, is negative. It is observed 
to be between and -0.3, whereas the model spectra have 
values in the range from -0.55 to -0.25 (Fig. [8j bottom) 
and thus are significantly more curved. 

In the cutoff distribution, the total fluxes of electrons 
escaping the acceleration site are initially 9.9 • 10 32 elec- 
trons s _1 above Et = 98 keV. This reduces to 6.6 • 10 31 
electrons s _1 for Et = 159 keV in the course of the de- 
cay. The total injected power reduces from 1.9 • 10 26 erg 
s _1 to 2.1 • 10 25 erg s _1 . In the turnover model, the total 
number of particles diminishes from 2.3 -10 33 to 1.3 -10 32 , 
and the injected power from 1.9 • 10 26 erg to 2.8 • 10 25 erg. 

6. DISCUSSION 

Spectroscopic RHESSI observations are well suited to 
study the spectral evolution during the main phase of 
large flares. The path observed in the 7 vs. log/50 plots 
for the events is not simple. However, it can be broken 
down reasonably into a superposition of linear sections 
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Fig. 7. — Photon spectra by thick target emission (at footpoints) 
from a constant power-law electron distribution having spectral 
index of <5 = —6.5 in the acceleration region. The low-energy cutoff, 
Et, increases from the top to the bottom curve from 0, 100, 120, 
140, to 160 keV. 



during flux rise and decay phases. While not all rise or 
decay phases can be so decomposed, this simple descrip- 
tion is adequate for most of them, and permits compari- 
son of observations and theory. 

There is a differen ce between the results reported here 
and the results from lGrigis fc Benz I (|2004f ) in the asym- 
metry between rise and decay phases in SHS peaks. The 
previous results indicated that decay phases are steeper 
in the 7 vs. log/50 plot than rise phases. We find the 
opposite. The reason is probably the selection bias: here 
we selected specifically events showing hardening. This 
hardening trend sometimes overlays SHS peaks, giving 
rise to a soft-hard-less-soft pattern. 

The hard X-ray images during the events show the 
usual morphology of hard X-ray solar flares: a low-energy 
coronal source and two or more high-energy footpoint 
sources. The position of the footpoint sources is strongly 
variable: it either moves smoothly or jumps from loca- 
tion to location. This reflects changes in the connection 
between the accelerator and the chromosphere, as well 
as in the location of the accelerator itself. 

The behavior observed in the images cannot be reduced 
to one simple scenario valid for all events. However, the 
observations suggest that there is no clear separation be- 
tween the SHS and the hardening phases: the former 
seems to smoothly merge into the latter. Even in the 
cases where the emission jumps at the onset of harden- 
ing (Table [2]), the footpoint behavior seems not to change 
radically. 

Alternatives to the scenario presented in Sect. [5] are 
conceivable. In particular, electron storage in the corona 
and slow release during the decay could be a possibility. 
As Coulomb interactions are faster at low particle en- 
ergies, the spectrum would harden with time while the 
released flux decreases. Noting that the hardening phase 
in the 19-JAN-2005 event lasts more than 30 minutes and 
that the decay of the flux in time is nearly exponential 
(as seen by the fact that the /50 line in Fig. [2j panel D, 
is nearly straight), the total number of injected electrons 
can be computed from the total electron fluxes at the 
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start -Fbeg an d a t the end -Fend- 

.Fend — Fbeg 



F 



TOT 



At- 



log 



END 



(4) 



F, 



BEG 



where At is the observed duration (here 30 minutes) and 
Ftot the total injected flux. From the observed values, 
we get Ftot — 1 • 10 36 electrons for the turnover model 
and Ftot — 6-10 35 electrons for the cutoff model. These 
numbers do not seem extraordinarily high, but it should 
be noted that all these electron have energies above 100 
keV. 

Would such a large population of electrons be seen as 
a coronal hard X-ray source in the 50 - 100 keV band? 
The luminosity depends on the volume and density in the 
storage region. The observed footpoints in the hardening 
phase are separated from each other by approximately 
60", indicating a medium sized loop. Therefore, it would 
be visible on RHESSI images unless it were excessively 
under-dense. Thus we reject the storage model. 

Both the cutoff and turnover model are able to repro- 
duce the observed 7 vs. log/50 trend, but fail to repro- 
duce the correct spectral curvature i] (Fig. [8]). Although 
the observation of r/ is more difficult in the decay phase 
due to the lower signal-to-background ratio, the differ- 
ence between the cutoff and turnover models and the 
observed points is significant. The value of the parame- 
ter r\ depends on the energy interval chosen for the fitting 
of the model photon spectra (20-80 keV in our case). A 
lower maximum energy of this interval produces lower 
model values for r\. 

We note furthermore that if the accelerator is inhomo- 
geneous, the electron spectrum at the footpoints is the 
superposition of different components with different val- 
ues of the low-energy cutoff or turnover Fx- The super- 
position of spectra that curve at different energies is less 
curved than individual spectra. Figure suggests that 
the superposition of components with different spectral 
shape may in fact reduce of the total spectral curvature. 
Thus we consider the disagreement in curvature not cru- 
cial to reject the model. 

7. CONCLUSIONS 

We have presented results of spectroscopy, imaging and 
simple modeling of spectral hardening observed occasion- 
ally in the hard X-ray emission of large flares. The main 
conclusions are: 

• The flares selected for the presence of a hard- 
ening phase also show soft-hard-soft behavior, at 
least initially. The hardening starts at or after the 
largest peak of the flares. In 3 out of 5 events it 
starts 2 to 6 minutes after the onset of a CME. 

• Similar to SHS peaks, hardening phases can usually 
be described by piecewise linear sections in a plot 
of spectral index vs. logarithmic flux. 

• There is no clear trend relating the behavior of hard 
X-ray footpoint sources with the spectral evolution 
that would be valid for all events. Sometimes the 
location of the emission shifts when the hardening 
starts, in other events it does not. 
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Fig. 8. — Observed values of 7 vs. /50 (top panel) and n vs. 
-Z50 (bottom panel), represented by asterisks. The best-fit model 
curves (continuous for the cutoff model, dashed for the turnover 
model) are expected from model with constant electron hardness 
and flux, but rising low energy cutoff or turnover energy Et- 

• In the event of 19-JAN-2005, there are only 
two well-defined footpoint sources throughout the 
whole event. No discontinuity is observed in the 
motion at the onset of hardening, but a general 
trend of slowing down, such that the FPs become 
nearly stationary during the decay phase. 

• In 3 out of 5 flares, the coronal source moved 
continuously during the onset of the hardening. 
This motion was directed upwards in two near-limb 
events. 

In the sample studied, we find a surprising lack of de- 
tailed correlation between spectral and s patial behavior. 
It is sim ilar to what has been observed bv lGrigis fc Benzl 
(2005b]) in a smaller flare (M6) featuring strong footpoint 
motions and hardening at the end. 

The main question addressed in this paper is whether 
the SHS peaks and the hardening phases are the results 
of two different acceleration mechanisms. The results 
support the view that the same acceleration mechanism 
changes gradually in the later phase of the flare. This 
change has clear effects on the spectrum, but a more in- 
direct and subtle influence on the source position. The 
operation of a second acceleration process later in the 
flare cannot be ruled out, however. Nevertheless, we have 
found strong evidence that there is a gradual change in 
the accelerator, transforming its behavior from impul- 
sive (showing up as SHS peaks) to gradual (hardening 
phases). This is substantiated by the observations of the 
superimposition of SHS peaks with a continuous hard- 
ening trend and of the smooth footpoint motions during 
the onset of hardening. 
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The reas on for the associa tion with interplanetary pro- 
ton events (Kiplingcr 1995) remains to be explored. As 
an aside we may note that when the footpoints drift 
apart, acceleration takes place in larger and larger loops. 
In a stochastic acceleration framework, the acceleration 
efficiency of electrons in larger loops is r educed, while 
ions can be more efficiently accelerated (jEmslie et al. I 
[Mil . Since hardening trends are well correlated with 
the occurrence of interplanetary energetic protons events, 
it is possible that the very conditions that are responsible 
for the hardening trends favor acceleration of protons, 
which may then escape into interplanetary space, with 
the CME controllin g their release rather than accelera- 
tion (jSimnett H20"oH . 

The observed motion of footpoints suggests that differ- 
ent coronal loops may be involved in particle acceleration 
during a flare. They will have different physical proper- 
ties such as size, density, and magnetic field. The overall 
magnetic geometry of the active region will determine 
which loops reconnect at which time, sometimes giving 
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rise to an orderly motion of footpoints, sometimes gener- 
ating a more chaotic situation. The data suggest that as 
the reconnection process proceeds, some physical param- 
eters of the acceleration site changes in such a way as to 
favor the production of harder spectra, rather than hav- 
ing a totally new process (say, shock acceleration) taking 
over in the decay phase. 
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APPENDIX 
HARD X-RAY MODEL FITTING 

The hard X-ray spectrum observed in solar flares consists of two distinct components at low energies (that is, below 
10-40 keV) and high energies. The properties of the low-energy component are typical of thermal emission of a hot 
plasma with temperature of 10 to 40 MK. In particular, the spectrum steepens with energy, and the temporal evolution 
follows the (cooler) thermal plasma observed in soft X-rays and EUV. 

The high-energy component behaves differently and is hence dubbed non-thermal. It is harder then the thermal 
component and is usually fitted with a power-law function of the energy with 2 free parameters. However, sometimes 
the observed spectru m steepens at higher energies. In the lit erature, this is usually accounted for by using a broken 
power- law model (e.g. lDulk et al. Ill992t iBattaglia et al. 1120051 ). There are some disadvantages in the broken power-law 
model: a) it is not physical, in the sense that any continuous electron distribution emitting X-rays should generate a 
differcntiablc photon spectrum, and b) the location of the break-point is poorly determined by the observations. 

We argue that there is a simpler extension to the power-law model which both turns down at higher energies and 
is smooth. Recalling that a power-law function plotted in log-log space is a straight line, we choose as a "natural" 
extension to the next order a parabolic model in log-log space, described by the function 

/ , \ 7 - V log (eAo) 
/(e) = I eo ■ -j . (Al) 

The 3 model parameters are the normalization I t01 amounting to the flux at the (fixed) normalization energy eo, the 
spectral index 7 < and the parabolic coefficient 7/, which we will refer to as the spectral curvature although, strictly 
speaking, the geometric curvature of the parabola is not constant, but equals — 2ri in the vertex and vanishes at infinity. 

In the special case 77 = 0, the (unbroken) power-law model is recover ed. We note here th at a log-parabolic model 
has been used previously to describe observed X-ray spectra of pulsars (Mas saro et al. Il2000f ) . 

In summary, the reasons for preferring the log-parabolic model over the more usual broken power-law are: 

1. It is simpler than the broken power-law, as it allows only 3 instead of 4 free parameters. 

2. For the vast majority of the time intervals, it produces similar values of x 2 as the broken-power law. 

3. It is diffcrentiable, therefore there exists a continuous electron spectrum producing the photon spectrum. This 
is not the case for the broken power-law, where a discontinuity is needed in the electron spectrum, which would 
quickly be eliminated by kinetic plasma processes. The spectral index increases linearly with loge. 

4. It allows a better comparison with acceleration models which naturally produce slightly curved electron spectra 
(like stochastic acceleration). 

Therefore, we fit the spectra to a photon model with an isothermal component at lower energies and a log-parabolic 
component as given above at higher energies. The background is taken into account in the following way: the pre- 
event and post-event background spectra are measured and averaged (in some cases, particle contamination prevented 
to obtain both of them, and only one was taken instead), yielding a reference background spectrum. The reference 
background spectrum is multiplied with a free parameter A and added to the model spectrum. A is fitted together 
with the other model parameters. 
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In large flares, som etimes an additional hard and weak 7-ray emission from electrons is observed above 200 keV 
(jKrucker et al J 12008( 1. In this paper we do not report on the properties of this emission since it is not related to the 
hardening of the spectrum at lower energies. Does this component negatively affects our fittings? Examination of the 
fitted spectra reveal that the fittings account for this weak, hard emission by an increase of the parameter describing 
the strength of the background (A) by a factor of about 1 to 2. The fittings are of good quality because the background 
is also weak and hard. This erroneous increase in the strength of the background has no effect at lower energies, since 
there the hard X-ray flux is stronger by orders of magnitude. 

PIVOT POINT AND PARABOLIC FITTING 

A linear dependence of 7 vs. log/ eo with negative slope can be interpreted geometrically as a fixed intersection point 
of the various power-law spectra at different times. The intersectio n in the spectral plane (/ vs. e) is called the pivot 
point, located below the reference energy eo (jGrigis fc Benz Il2005all . Similarly, a linear dependence with positive slope 
can be interpreted by a pivot point at energy larger than eo- If fitting a log-parabolic spectrum, it is the tangents to 
the spectrum at eo m log-log space that are intersecting, rather than the curves themselves. 

The description in terms of a pivot point has the advantage that it does not depend on the choice of the reference 
energy eo- O n the other hand, the pivot point energy jumps from to +00 when the slope in 7-log/ eo goes from 
negative to positive. The relation between the pivot point coordinates e»,/* and the line parameters a, b in Eq. [T]are 
given by 

a = 1 b= -log/* , m , 

log(e*/e ) log(e»/e ) ' 

Since the spectra are curved in log- log space, the local spectral index 7', defined as the logarithmic derivative of the 
spectrum 

= g log /(e) = _ ^ 
aloge 

is energy dependent. The spectral parameter 7 is equal to the local spectral index 7' at the reference energy eo = 50 keV. 
The presence of a strong correlation of the time series of 7 and log /50 does not necessarily imply a strong correlation 
of the time series of Y(e) and log /(e) at energies e ^ e . In fact, if the spectra have a common pivot point at e* < e , 
the correlation of 7'(e) and log 1(e) weakens near the pivot point, and turns into an anticorrelation for e < e*. On 
the other hand, if a pivot point exist at e* > e , then the anticorrelation between 7'(e) and log /(e) for e ~ e < e* 
transforms itself into a correlation at e > e*. 

An examination of the data for the events studied here shows that in the time intervals when the correlation 
coefficient between log/50 and 7 is near 1, the correlation coefficient between 7'(e) and log /(e) is approximately 
constant at energies higher than 30-40 keV However, it shows a strong decrease to values close to -1 at lower energies 
such that the transition takes place around 5-20 keV, near the pivot point position. On the other hand, when the 
correlation coefficient is near -1 at 50 keV we observe that the correlation coefficients rises toward 1 at higher energies, 
as these events tend to have a pivot point at energies larger than 50 keV. 

Therefore, the energy dependence of the correlation in the data follows a similar pattern as the one expected for 
non-curved spectra. This is because the observed spectra are not strongly curved: rj < 0.257 in all the fitted spectra, 
and 77 < 0.157 in 85% of all fitted spectra. 

COMPARISON BETWEEN MODEL AND OBSERVATION 

Comparison of non-thermal hard X-ray spectral observations and theoretical models can be performed in different 
ways. In our case, we have a really simple model and a time-dependent situation. The goal of the comparison is not 
the perfect reproduction of every single observed spectrum, but rather a coherent description of the time evolution 
which should be compatible with the observed data. This is implemented by the additional constraint in the model (as 
given by Eq. [2]) of keeping a constant electron spectral index 5 above the threshold energy Et- So our task consists of 
two steps: first, for a given value of 8, we have to find the set of E T (t) and F Eo (t) that best reproduces the data and 
second, we have to choose the best value of S (for instance by running the first step for many different values of S and 
pick the best one). 

For the comparison of the model with the data, we generate model photon spectra emitted by the model electron 
spectra by computing the thick-target Bremsstrahlung emissi on assumi ng collis i onal e nergy losses and using the full 
relativistic Bethe-Heitler cross section (|Bethe fc Heitler"l ll934) with the lElwertl (|1939f ) correction factor. 

For the first step, the obvious strategy involves forward fitting the electron model to the observed data. However, 
this is time consuming because it needs to be repeated for all the spectra and many different values of 6. On the other 
hand, we have at our disposal the photon fitting parameters / eo , 7 and 77 (Eq. IA1|) . which are good descriptions of the 
observed photon spectrum as long as the reduced x 2 is around one. In this case, there is no need of additional fitting in 
count space: we just fit exactly the same log-parabolic model to the model photon spectra (in the energy range where 
the non-thermal component is seen above the thermal component and the background in the observations). This is 
faster, and delivers photon model parameters I™ , -y MOD and 7y MOD as a function of the electron model parameter 
Fe i 8 and Et- The comparison can then be performed in the 7 vs. I eo plot by a least square argument in two steps. 

In the first step we held 5 fixed and increase Et to generate a curve in the 7 vs. / eo plot for each value of Fe - The 
normalization Fe is then chosen such that it minimizes the square differences between / £o and I^ OB - The results 
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from running step one repeatedly are a set of paths (one for each different value of S) in the 7 vs. I fQ plot, where the 
variable E T changes along the path. 

In the second step, we just select among all paths in the 7 vs. I eo space the one with the least square distances from 
all the observed points. This yields the best 5 and Fe and a range of variation of Et- 
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